close all

load fig_number

    
IRF_periods = 20;


if figure4==1
load IRF_baseline

fig1 = figure;
subplot(3,3,1)
plot(Y_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Output','FontSize',14)
grid on

subplot(3,3,2)
plot(C_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Total Consumption','FontSize',14)
grid on

subplot(3,3,3)
plot(C_m_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Market Consumption','FontSize',14)
grid on

subplot(3,3,4)
plot(N_m_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Market Hours','FontSize',14)
grid on

subplot(3,3,5)
plot(N_n_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Non-Market Hours','FontSize',14)
grid on

subplot(3,3,6)
plot(W_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Wage','FontSize',14)
grid on


subplot(3,3,7)
plot(MU_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Markup','FontSize',14)
grid on


subplot(3,3,8)
plot(Pi_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Inflation','FontSize',14)
grid on


if alpha1I==0
   subplot(3,3,9)
   plot(SigmaA_net(1:IRF_periods),'b','LineWidth',1.5); hold on
   plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
   title('Exogenous Demand Volatility','FontSize',14)
   grid on
   fig_name = 'IRF_Baseline';   

elseif alpha1I==1 

   subplot(3,3,9)
   plot(Alpha1_net(1:IRF_periods),'b','LineWidth',1.5); hold on
   plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
   title('Market Consumption Share','FontSize',14)
   grid on
   fig_name = 'IRF_Alpha';   
     
end
    
end

if figure3==1

load IRF_baseline IRF_no_homeprod

fig1 = figure;
subplot(3,3,1)
plot(Y_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(Y_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Output','FontSize',14)
grid on

subplot(3,3,2)
plot(C_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(C_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Total Consumption','FontSize',14)
grid on

subplot(3,3,3)
plot(C_m_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(C_m_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Market Consumption','FontSize',14)
grid on

subplot(3,3,4)
plot(N_m_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(N_m_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Market Hours','FontSize',14)
grid on

subplot(3,3,5)
plot(N_n_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(N_n_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Non-Market Hours','FontSize',14)
grid on

subplot(3,3,6)
plot(W_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(W_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Wage','FontSize',14)
grid on


subplot(3,3,7)
plot(MU_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(MU_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Markup','FontSize',14)
grid on


subplot(3,3,8)
plot(Pi_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(Pi_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Inflation','FontSize',14)
grid on


subplot(3,3,9)
plot(SigmaA_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(SigmaA_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Exogenous Demand Volatility','FontSize',14)
grid on
fig_name = 'IRF_Baseline';   

legend('Home Production','No Home Production')
legend('boxoff')  

end

if figure5==1 
IRF_periods = 15;
    
fig1 = figure;
subplot(3,2,1)
plot(Y_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(Y_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Output','FontSize',14)
grid on

subplot(3,2,2)
plot(C_m_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(C_m_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Market Consumption','FontSize',14)
grid on


subplot(3,2,3)
plot(N_m_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(N_m_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Market Hours','FontSize',14)
grid on


subplot(3,2,4)
plot(N_n_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(N_n_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Non-Market Hours','FontSize',14)
grid on


subplot(3,2,5)
plot(W_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(W_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Wage','FontSize',14)
grid on



subplot(3,2,6)
plot(Pi_net(1:IRF_periods),'b','LineWidth',1.5); hold on
plot(Pi_net_nohp(1:IRF_periods),'r--','LineWidth',1.5); hold on
plot(zeros(IRF_periods,1),'k--','HandleVisibility','off'); xlim([1 IRF_periods]);set(gca,'XTick',[4:4:IRF_periods],'FontSize',12);
title('Inflation','FontSize',14)
grid on

legend('Uncertainty (No home production)','Uncertainty and Reallocation Policy','Location','SouthEast')
legend('boxoff')  

fig_name = 'IRF_Covid';   
    
end



% cd D:\Dropbox\Reserach\Papers\Uncertainty_and_Housework\Text\Paper\Figures
% set(gcf, 'Position', get(0, 'Screensize'));
% orient(fig1,'landscape')   
% print(fig1, '-dpdf', [fig_name])
% 
